ANHARMONIC OSCILLATORS IN THE COMPLEX PLANE, 
PT-SYMMETRY, AND REAL EIGENVALUES. 
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Abstract. For integers m > 3 and 1 < £ < m — 1, we study the eigenvalue problems 
—u"(z) + [(— l) l (iz) m — P{iz)]u(z) = Au(z) with the boundary conditions that u(z) decays 
to zero as z tends to infinity along the rays argz = — -| ± m the complex plane, 

where P is a polynomial of degree at most m — 1. We provide asymptotic expansions of the 
eigenvalues A„. Then we show that if the eigenvalue problem is PT-symmetric, then the 
(3J^' eigenvalues are all real and positive with at most finitely many exceptions. Moreover, we 

show that when gcd(m,£) — 1, the eigenvalue problem has infinitely many real eigenvalues 
if and only if its translation or itself is 'P7"-symmetric. Also, we will prove some other 
, interesting direct and inverse spectral results. 

2010 Mathematics Subject Classification: 34L4O, 34L20, 81Q12 
Key words: Anharmonic oscillators, asymptotics of the eigenvalues, VT -symmetry 

-1— > 

a . 

1. Introduction 

In this paper, we study Schrodinger eigenvalue problems with real and complex polynomial 

>• ' potentials in the complex plane under various decaying boundary conditions. We provide 
IT) ' 

O I explicit asymptotic formulas relating the index n to a series of fractional powers of the 
eigenvalue X n (see Theorem 1.1). Also, we recover the polynomial potentials from asymptotic 
formula of the eigenvalues (see Theorem 1.4 and Corollary 2.3) as well as applications to the 
§ so-called PT-symmetric Hamiltonians (see Theorems 1.2 and 1.3). 

For integers m > 3 and 1 < £ < m — 1, we consider the Schrodinger eigenvalue problem 

d 2 



X : (i-i) (H iU 



+ (-iy(iz) m - p(iz) 



u{z) = Xu(z), for some A G C, 



dz 2 

with the boundary condition that 

7T (£ + lW 

(1.2) u(z) — > as z — > 00 along the two rays argz = ± , 

K J K J 6 J 6 2 m + 2 ' 

where P is a polynomial of degree at most m — 1 of the form 

P(z) = a 1 z m ~ 1 + a 2 z m ~ 2 H h a m _iz + a m , ^ G C for 1 < j < m. 

If a nonconstant function u satisfies (1.1) with some A G C and the boundary condition 
(1.2), then we call A an eigenvalue of He and u an eigenfunction of Hi associated with the 
eigenvalue A. Sibuya [27] showed that the eigenvalues of He are the zeros of an entire function 
of order p := 0+ (0) 1) an d hence, by the Hadamard factorization theorem (see, e.g., [7, p. 
291]), there are infinitely many eigenvalues. We call the entire function the Stokes multiplier 



Date: August 6, 2010. 



2 

(or the spectral determinant), and the algebraic multiplicity of an eigenvalue A is the order 
of the zero A of the Stokes multiplier. Also, the geometric multiplicity of an eigenvalue A is 
the number of linearly independent eigenfunctions associated with the eigenvalue A, that is 
1 for every eigenvalue A [16, §7.4]. 

We number the eigenvalues {A n } n >jVo in the order of nondecreasing magnitudes, counting 
their algebraic multiplicities. We will show that the magnitude of large eigenvalues is strictly 
increasing (see, Lemma 2.2) and hence, there is a unique way of ordering large eigenvalues, 
but this is not guaranteed for small eigenvalues. However, how we order these small eigen- 
values will not affect results in this paper. Throughout this paper, we will use A n to denote 
the eigenvalues of Hg without explicitly indicating their dependence on the potential and the 
boundary condition. Also, we let 

a : = (ai,a 2 , . . . , a m ) G C m 

be the coefficient vector of P{z). 

The anharmonic oscillators Hg with the various boundary conditions (1.2) are considered 
in [4, 24]. When m is even and i — y, H™ is a Schrodinger operator in L 2 (R) (see, e.g., 
[1, 2, 5, 11, 12, 15, 18, 19]). This is self-adjoint if the potential V{z) = (-lY(iz) rn - P(iz) 
is real on the real line, and non-self- adjoint if the potential is non-real. 

Some particular classes of Hi have been studied extensively in recent years in the context 
of theory of PT-symmetry [3, 6, 10, 14, 20, 26, 28, 29]. The Hg is PT-symmetric if the 
potential V satisfies V(—z) = V(z), z G C, that is equivalent to a G W rn . In this paper, we 
will generalize results in [26] (where Hi is studied) tol<^<m — 1 and introduce some 
new results. These results are consequences of the following asymptotic expansion of the 
eigenvalues. 

Theorem 1.1. For each integer m > 3 and 1 < i < m — 1, there exists an integer Nq = 
N (m,t) such that the eigenvalues {A„}„>at of He satisfy 

1 m+1 i h± 

(1.3) n + - = J2 c ej( a ) X n +m+r le(a) + 0(\- p ) as n ->■ oo, 

3=0 

where Qj(a) and rje(a) are defined in (2.4) and (2.5), respectively. 

Also, we obtain the partial reality of the eigenvalues for PT-symmetric Hg. 

Theorem 1.2. If Hi is VT -symmetric, then eigenvalues are all real with at most finitely 
many exceptions. 

Proof. If u(z) is an eigenfunction associated with the eigenvalue A of a "PT-symmetric Hg, 
then u(—z) is also an eigenfunction associated with A. In Corollary 2.2, we will show that 
I A n | < |A n +i| for all large n and arg(A n ) — > 0. Thus, A n are real and positive for all large n 
since |A| = |A|. □ 
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Many "PT-symmetric operators have real eigenvalues only [10, 22, 23]. However, there are 
some PT-symmetric He that produce a finite number of non-real eigenvalues [2, 8, 9, 14]. 
When Hi is self-adjoint, the spectrum is real. Conversely, when the spectrum is real, what 
can we conclude about Hp. The next theorem provides a necessary and sufficient condition 
for Hi to have infinitely many real eigenvalues. 

Theorem 1.3. Suppose that gcd(m, f) = 1. Then H e with the potential V(z) = — (iz) m + 
P{iz) has infinitely many real eigenvalues if and only if Hi with the potential V(z — z ) for 
some zq G C is VT ' -symmetric. 

The next theorem reveals an interesting feature of the eigenvalues as a sequence: 

Theorem 1.4. Suppose that gcd(m, £) = 1. Let {A n }„>jv and {A n }„>Ar be the eigenvalues 
of Hi with the potentials V and V, respectively. Suppose that X n — X n = o(l) as n — >■ oo. 
Then V{z) = V(z — z ) for some z G C and X n = X n for all n > N after, if needed, small 
eigenvalues are reordered. 

The asymptotic expansions of the eigenvalues of Hi with £ = [ttJ have been studied in, for 
example, [13, 15, 18]. Maslov [18] computed the first three terms of asymptotic expansions 

3 2 

of Xn, where X n are the eigenvalues of — J^w + x A u = Xu, u G L 2 (R). Helffer and Robert 
[15] considered 

d 2k 

--^u + (x 2m +p(x))u = Xu, ueL 2 (R), 

where k, m are positive integers and where p(-) is a real polynomial of degree at most 2m — 1. 
They obtained the existence of asymptotic expansions of the eigenvalues to all orders, and 
suggested an explicit way of computing the coefficients of the asymptotic expansion. In 
particular, for the case when the potential is ex^+x 2 , e > 0, Helffer and Robert [15] computed 
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the first nine terms of the asymptotic expansion of A^. Also, Fedoryuk [13, §3.3] considered 
(1.1) with complex polynomial potentials and with (1.2) for £ = [yj and computed the 
first term in the asymptotic expansion. Also, Sibuya [27] computed the first term in the 
asymptotic expansion for £ = 1. 

This paper is organized as follows. In Section 2, we define cgj(a) and some other notations. 
Also, we invert (1.3), expressing A n as a series of fractional powers of the index n and prove 
Theorems 1.3 and 1.4, and other interesting direct and inverse spectral results. In Section 
3, we introduce some properties of the solutions of the differential equation in (1.1), due 
to Hille [16] and Sibuya [27]. We study the asymptotic of the Stokes multiplier associated 
with Hi in Section 4 and treat the general case H e in Section 5. In Section 6, we relate 
the eigenvalues of He with the zeros of the Stokes multiplier. We prove Theorem 1.1 for 
1 < £ < f in Section 7 and for £ > f in Section 8. 



2. Notations and corollaries 



In this section, we will define cgj(a) and some other notations and introduce some corol- 
laries of Theorem 1.1. 

We define, for nonnegative integers k,j, 



(2.1) 



bjJa) is the coefficient of 



M8( 



a>i a 2 a m \ k 
1 — 7T -\ 1 and 

z z 2 z m J 



(2.2) 



fc=0 



Notice that &o,o( a ) — &o(°) = 1, bj >0 (a) — if j > 1, and bj^{a) = if j < k or k < ^. We 
will also use 



(2.3) 



K m ,j{a) = ^2 K m,i,k bj,k(a) for j > 0, 



fc=0 



where for k < i-, K m ,- fc = and for ^- < < ?, 



Km,j,k ■ < 



2 cos -2L 

V m 
__2_ 

m 



if j = fc = 0, 

if J = fc = l, 



_2_ 

m 



m even, 



(ln2-Eti 1 2 l L T) ifi = f + l, 
I ifjVlorjVf + 1- 



Now we are ready to define Qj(a) as follows: for 1 < j < m + 1, 

j 



k=0 



(2.4) Q » = — £(-l) ( ' +1)fc /W*W«) 
and we also define 

(2.5) 7fc(a) 
where 

i/(a) 



1 cos 



m 



_ 1) ^Ma) if£isodd 

/ m 



if otherwise 



m 



and fi(a) — — — v(a), 



bm. + i(a) if m is even, 
if m is odd. 



2.1. Further direct spectral results. One can invert the asymptotic formulas (1.3) to 



obtain formulas for A n in terms of n. 



Corollary 2.1. One can compute dej(a) explicitly such that 

(2.6) A n = $><j(a)U + -J +0(n-^j, 

_ 2m 

where d efi (a) = (n^B (1, 1 + sin (|)) ™+ 2 > 0. 

Proof. Equations (1.3) is an asymptotic equation and it can be solved for A n , resulting in 
(2.6). For details, see, for example, [24, §5]. □ 

In the next corollary, we provide an asymptotic formula for the nearest neighbor spacing 
of the eigenvalues. And the large eigenvalues increase monotonically in magnitude and have 
their argument approaching zero. 

Corollary 2.2. The space between successive eigenvalues is 

( m-2 
1\ m + 2 / m-2 \ 

n+- +o . 
2/ V / 

In particular, lim^^ |A n+ i — A n | = oo and lim n _>. +00 arg(A n ) = 0. Hence: 

I A n | < |A n+ i| for all large n. 

Proof. These claims are consequences of (2.6) and the generalized binomial expansion. For 
details, see, for example, [24, §6]. □ 

Suppose that X n are eigenvalues of Hi for some P. Then the degree m of the polynomial 
potential can be recovered by 

2m _ lim l^ <2 ^ 



m + 2 n-5>oo In n 

The next corollary shows that one can recover the polynomial potential from the asymptotic 
formula for the eigenvalues. 

Corollary 2.3. Suppose that gcd(m,£) = 1 and that {A„}„>jv are the eigenvalues of H e 
with the potential V . If 

m+l . / 1 \ 

(2.8) C *A + ^ + O (a^~™) = In + - J 
as n — y oo for some c* G C, or if 

(2.9) X n = J2d*(n+-\ +0(n-—iy 

3=0 ^ ' 

for some d* G C 7 then there exists V such that for each polynomial V with which Hi generates 
{A n }„>jv , V(z) = V (z - z ) for some z G C. 

Proof. The coefficients Qj(a) and dgj(a) have the following properties. 
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(i) The cej(a) and d^j{a) are all real polynomials in terms of the coefficients a of 
P{x). 

(ii) The coefficients ce t o(a) and de t o(a) do not depend on a (they are constants). 

(iii) For 1 < j < m, the polynomials cij(a) and dej(a) depend only on a±, 0,2, ■ ■ ■ , 
Furthermore, if gcd(m, £) = 1 and if j 7^ 1, y + 1, then Qj(a) and dej(a) are non- 
constant linear functions of a,-. 

Here we will sketch the proof when (2.8) is assumed. First, notice that c\ = 0. Otherwise, 
\ n are not the eigenvalues of Hg, according to (2.4). Choose a\ = 0. Then since Q^a) is a 
non-constant linear function of a<i and depends on a\ and 02, 0^2(0) = C2 determines a-i- If m 
is even, Q > m +1 (a) = and we use 7^(a) instead of c^+^a). Suppose that for 2 < j < m — 1, 
the C2, . . . , c* uniquely determine 02, 03, . . . ,Oj. Then by (iii), Q J+ i(a) = c* +1 determines 
Oj+i uniquely. So by induction, one can determine all cij for 2 < j < m and Vo is the potential 
with these a/s. Then V(z) = Vo(z — zq) for some zq G C. Otherwise, the eigenvalues {A„} 
do not satisfy (2.8). 

The case when (2.9) holds can be handled similarly. □ 

Next, we will prove Theorems 1.3 and 1.4. 

Proof of Theorems 1.3. If V(z — z ) for some z Q e C is PT-symmetric, then by Theorem 1.2, 
all but finitely many eigenvalues are real and hence, there are infinitely many real eigenvalues. 

Suppose that with the potential V(z) has infinitely many real eigenvalues. Then we 
can always find z G C so that V(z — z ) has no ^ m ~ 1 -term, that is, a± = 0. We will show 
that V(z — z ) is PT-symmetric. Since there are infinitely many real eigenvalues, C£j(a), 
< j < m + 1, and r]e in (1.3) are real. Also, the sine term in (2.4) does not vanish except 
for j = 1. And the cosine term does not vanish except for j = y + 1 when m is even in 
which r]e replaces the role of c^m. 

Like we did for proof of Corollary 2.3, since a\ = 0, from the above properties of Qj and 
by the induction, we can show that aj for 2 < j < m are all real and hence, V(z — z ) is 
"PT-symmetric. □ 

Proof of Theorems 1.4. Suppose that X n — \ n = o(l) as n — > 00. Then dij(a) = di^ia) for 
all 1 < j < ra + 1. Then Corollary 2.3 completes proof. □ 

3. Properties of the solutions 

In this section, we introduce work of Hille [16] and Sibuya [27] about properties of the 
solutions of (1.1). 

First, we scale equation (1.1) because many facts that we need later are stated for the 
scaled equation. Let u be a solution of (1.1) and let v(z) = u(—iz). Then v solves 

(3.1) -v"{z) + [(-lY +1 z m + P(z) + \}v(z) = 0. 
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Figure 1. The Stokes sectors for m = 3. The dashed rays represent arg^ = ±|, i^f , 7r. 
When £ is odd, (3.1) becomes 

(3.2) -v"(z) + [z m + P(z) + A]u(z) = 0. 

Later we will handle the case when £ is even. 

Since we scaled the argument of u, we must rotate the boundary conditions. We state 
them in a more general context by using the following definition. 

Definition. The Stokes sectors S k of the equation (3.2) are 

2kn 



S k = UeC: 



arg{z) 



< — ^— \ for k E Z. 

m + 2 



m + 2 
See Figure 1. 

It is known from Hille [16, §7.4] that every nonconstant solution of (3.2) either decays to 
zero or blows up exponentially, in each Stokes sector Sk- 

Lemma 3.1 ([16, §7.4]). 

(i) For each k G Z, every solution v of (3.2) is asymptotic to 



(3.3) (const. )z 4 exp 



as z — > oo in every closed subsector of S k . 
(ii) If a nonconstant solution v of (3.2) decays in Sk, it must blow up in Sk-i U Sk+i- 
However, when v blows up in Sk, v need not be decaying in Sk-i or in Sk+i- 

Lemma 3.1 (i) implies that if v decays along one ray in Sk, then it decays along all rays in 
Sk- Also, if v blows up along one ray in Sk, then it blows up along all rays in Sk- Thus, the 
boundary conditions (1.2) with 1 < £ < m — 1 represent all decaying boundary conditions. 

Still with £ odd, the two rays in (1.2) map, by z (->■ — iz, to the rays arg(z) = ± ^qjrr 
which are the center rays of the Stokes sectors Se+i and S i+i and the boundary conditions 



(1.2) on u become 

v decays to zero in the Stokes sector Se+i and S e+i . 

2 2 

When £ is even, we let y{z) = v{u~^z) so that (3.1) becomes 

(3.4) -v"(z) + [z m + u^Piuj-h) + W 1 X]y{z) = 0, 
where 

2m 

u = exp 

F [m + 2_ 

and hence, co~^ +1 = —1. For these cases, the boundary conditions (1.2) become 
y decays to zero in the Stokes sector Se+2 and S_t. 

2 2 

The following theorem is a special case of Theorems 6.1, 7.2, 19.1 and 20.1 of Sibuya [27] 
that is the main ingredient of the proofs of the main results in this paper. For this we will 
use r m = — 22 if m is odd, and r m = — ^ — &™ +1 (a) if m is even. 

Theorem 3.2. Equation (3.2), with a G C m , admits a solution f(z,a,X) with the following 
properties. 

(i) f(z, a, A) is an entire function of z, a and A. 

(ii) f(z, a, A) and f'(z, a, A) = -§^f(z, a, A) admit the following asymptotic expansions. 
Let e > 0. Then 

f(z, a, A) = z r -(l + Oiz- 1 ' 2 )) exp [-F(z, a, A)] , 
f(z, a, A) = - + 0(^ 2 )) exp [-F(z, a, A)] , 

as z tends to infinity in the sector \ aigz\ < — e, uniformly on each compact set 
of (a, X) -values . Here 

F(z,a,X) = -^—z^ +1 + Y 1 bj(a)z^ m+2 - 2j \ 

(iii) Properties (i) and (ii) uniquely determine the solution f(z,a,X) of (3.2). 

(iv) For each fixed a G C m and S > 0, f and f also admit the asymptotic expansions, 

(3.5) /(0, a, A) =[1 + O (X- p )]X- lfi exp [L(a, A)] , 

(3.6) /'(0, a, A) = - [1 + O (A-^A 1 / 4 exp [L(a, A)] , 

as X ^ oo in the sector | arg(A)| < tx — 5, uniformly on each compact set of a G C m , 
where 

j Jo °° (\A m + p (t) + A - *t - £Tjf 6j(a)tT-^ i/m odd, 
L(a ' A) ~ j J +o ° (v/t™ + P(t) + A - tT - YJ =1 bj(a)tf-i - "-^f) dt if m is even. 

(v) The entire functions X >->■ /(0, a, A) and A >->■ /'(0, a, A) nave orders \ + —■ 
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Proof. In Sibuya's book [27], see Theorem 6.1 for a proof of (i) and (ii); Theorem 7.2 for a 
proof of (iii) with the error terms o(l); and Theorem 19.1 for a proof of (iv). Moreover, (v) is 
a consequence of (iv) along with Theorem 20.1 in [27]. The error terms in (iii) are improved 
from o(l) to O (A~ p ) in [25]. Note that properties (i), (ii) and (iv) are summarized on pages 
112-113 of Sibuya [27]. □ 

Remarks. Throughout this paper, we will deal with numbers like (u a \) s for some s G K 
and «6C. As usual, we will use 

2ni 



u a = exp 



a- 



m + 2 



27T 

m + 2 



s g 



and if arg(A) is specified, then 

arg ((co a X) s ) = s [arg(w a ) + arg(A)] = s 
Lemma 3.3. Let m > 3 and a G C m be fixed. Then 

oo , - 

(3.7) L(a, A) = £ K m »A^ - ^ ln(A) 

as A — >■ oo in the sector \ arg(A)| < ir — 5, uniformly on each compact set of a e C m . 

Proof. See [24] for a proof. □ 

Sibuya [27] proved the following corollary, directly from Theorem 3.2, that will be used 
later in Sections 4 and 5. 

Corollary 3.4. Let a G C m be fixed. Then L(a, A) = K m \^ + ™{1 + o(l)) as A tends to 
infinity in the sector \ argA| < 7T — 5, and hence 

(3.8) Re (L(a, A)) = K m cos arg(A)) |A|^(1 + (1)) 

as X oo in the sector | arg(A)| < 7r — 5. 

In particular, Re (L(a, A)) — >■ +oo as A — >■ oo in any closed subsector of the sector 

|arg(A)| < -^2- In addition, Re (L(a, A)) — > — oo as A — > oo m any closed subsector of 
the sectors < I arg(A) | < ty - 5. 

Based on Corollary 3.4, Sibuya [27, Thm. 29.1] also computed the leading term in (1.3) 
for t — 1. Also, Sibuya [27] constructed solutions of (3.2) that decays in Sk, k G Z. Before 
we introduce this, we let 

(3.9) G e (a) := (u (m+1)e a u u me a 2 , • • • , u 2i a m ) for t G \l. 
Then we have the following lemma, regarding some properties of G e (-). 
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Lemma 3.5. For a G C m fixed, and l u t 2 ,le |Z, G tl (G l *(a)) = G h+e2 (a), and 

b hk {G l {a)) = cu« m+ V k -i%, k (a), i G ±Z. 

In particular, 

bj{G\a)) =u~%{a), I G Z, 

Next, recall that the function f(z, a, A) in Theorem 3.2 solves (3.2) and decays to zero 
exponentially as z — > oo in So, and blows up in S-i U Si. One can check that the function 

f k (z,a,\) := f(co- k z,G k (a),co 2k \), k G Z, 

which is obtained by scaling f(z,G k (a),u 2k X) in the ^-variable, also solves (3.2). It is clear 
that fo(z,a,\) = f(z,a,X), and that f k {z,a,X) decays in S k and blows up in S k -i U S k+ i 
since f(z, G k (a),u 2k X) decays in 5 - Since no nonconstant solution decays in two consecutive 
Stokes sectors (see Lemma 3.1 (ii)), f k and f k+ i are linearly independent and hence any 
solution of (3.2) can be expressed as a linear combination of these two. Especially, for each 
k E Z there exist some coefficients C k (a, A) and C k (a, A) such that 

(3.10) f k (z, a, A) = C k (a, X)f (z, a, A) + C k (a, X)f-i(z, a, A). 
We then see that 

/on\ n i \\ Wfe-i(a,A) ~ W fci0 (a,A) 

(3.11) C k (a, A) = — — — ■ — - — -— and C k (a, A) - 



W_i,o(a,A) KV ' 7 W_i )0 (a,A)' 

where VVj^ = fjfg — fjfe is the Wronskian of fj and Since both fj, fi are solutions of the 
same linear equation (3.2), we know that the Wronskians are constant functions of z. Also, 
fk and f k+ i are linearly independent, and hence W k , k +i ^ for all k G Z. 

Also, the following is an easy consequence of (3.10) and (3.11). For each k, £ G Z we have 

VMa, A) = C k (a, \)W t>0 {a, A) + C k (a, X)W t ,- 1 (a, A) 
, , Wfc _i(a, A)Wi i0 (a, A) W M (a, A)W € ,-i(a, A) 

1 J " W_i >0 (a,A) + W_i >0 (a,A) 

Moreover, we have the following lemma that is useful later on. 
Lemma 3.6. Suppose k, j G Z. Then 

(3.13) Wk+i,j+i(a, A) = w-'Wfcj^ajVA), 

and W ,i(a, A) = 2u^ a) , where fi(a) = f - i/(a). 



Proo/. See Sibuya [27, pp. 116-118]. 



□ 
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4. Asymptotics OF W_i,i(a, A) 

In this section, we introduce asymptotic expansions of W_i i i(a, A) as A — > oo along the 
rays in the complex plane [24]. 

First, we provide an asymptotic expansion of the Wronskian Woj(a,A) of f and fj that 
will be frequently used later. 

Lemma 4.1. Suppose that 1 < j < y + 1. Tnen for each a G C m ; 

(4.1) W 0J (a, A) = [2iw-4 + O (A - ')] exp [L(G J (a), u 2j ~ m - 2 X) + L(a, A)] , 
as A — >■ oo in the sector 

(4.2) -7T + 5 < 7T - — ] — + 5< arg(A) < tt - o". 
v ' ~ m + 2 ~ 

Next, we provide an asymptotic expansion of W_i ; i(a, A) as A — > oo in the sector near the 

negative real axis. 

Theorem 4.2. Let m > 3, a G C m and < 5 < ^ be fi xed - Then 

(4.3) W_i,i(o, A) = [2i + O (\- p )] exp [L(G _1 (a), w _2 A) + L(G(a), uj- m X)} , 
as A — > oo a/ong £ne rays in the sector 

(4.4) 7i h 5 < arg(A) < tt H 5. 

v 7 m + 2 - v y - m + 2 

Proof. This is an easy consequence of Lemma 4.1 with j = 2 and (3.13). □ 

Also, for integers m > 4 we provide an asymptotic expansion of W_i 5 i(a, A) as A — > oo in 
the sector | arg(A) | < ix — S. 

Theorem 4.3. Let a G C m and < 5 < 2(r7 f +2) be fixed. If m > 4 then 

W_i,i(a,A) =[2w5+^) + 0(A-")]exp [L(G _1 (a), uT 2 A) - L(a, A)] 

(4.5) + [2a;5+^)+ 2 ^) + (A^)] exp [L(G(a),u 2 X) - L(a, A)] , 
as A — )> oo in the sector 

(4.6) -tt + S < arg(A) < tt - 5. 

Next, we provide an asymptotic expansion of W-i,i(a, A) as A — > oo along the rays in the 
upper- and lower- half planes. 

Corollary 4.4. Let m > 4, a G C m and < 5 < ^ &e fi xed - Then 

W_i,i(a,A) = [2a;5+M«) +0(A^)]exp [L(G _1 (a), w _2 A) - L(a, A)] , 
as A — >■ oo in the sector 5 < arg(A) < n — 5. Also, 

W-!,i(a, A) = [2o;3+p(«)+M«) + ( A -/>)] exp [L(G(a), u 2 X) - L(a, A)] , 
as A — >■ oo in the sector —n + 5 < arg(A) < —5. 
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Proof. We will determine which term in (4.5) dominates in the upper and lower half planes. 
Since, by (3.8), 

Re (L(a, A)) = K m cos (^±^arg(A)) |A|^(1 + o(l)), 



we have 



[Re(L(G~ 1 (a),uj~ 2 X)) — Re(L(a, A))] - [Re (L(G(a), uj 2 X)) - Re (L(a, A))] 



27r m + 2 
cos ( + — 

m 2m 



{2n m + 2 ,,,1,1,, /1N , 

arg(A)J - cos ^— + arg(A) J j |A|* + ~(l + o(l)) 



= 2K m sin sin (^arg(A)) |A|^(l + o(l)). 

Thus, the first term in (4.5) dominates as A — > oo along the rays in the upper half plane, 
and the second term dominates in the lower half plane. This completes the proof. □ 

Proof of Theorem 4-3. In [24], C(a, A) is used for ^j^^y and asymptotics of C(a, A) are 
provided. Notice that W_i,i(a, A) = 2w^ a )C(a, A). 
Theorem 13 in [24] implies (4.5) for the sector 

41 -|7T 47T 

(4.7) ^-^- + 5<arg(A)<7r-— — -6. 

m + 2 m+2 

Theorem 14 in [24] implies that 

W_i,i(a, A) =[2u^ a) + O (\- p )] exp [L(G _1 (a), uT 2 A) - L(o, A)] 

as A — > oo in the sector it — +^ + 5 < arg(A) < n — 5. One can check that the first term 
dominates in this sector, by using an argument similar to that in the proof of Corollary 4.4. 
Also, Theorem 15 in [24] implies that 

W-i,i(a, A) =[2u 1+2 ^ a) + O (A"')] exp [-L(a, uj- m - 2 X) - L(- 2 (a), uT 4 A)] 

+ [ 2a; 3+M«)+M«) + o (x~P)] e X p [L(G(a),uT m A) - L(a,uj- m ~ 2 \)] , 

as A — > oo in the sector it + 5 < arg(A) < n + ~~ O ne can check that the second 
term dominates in this sector. Then we replace A by u m+2 X to convert the sector here to 
—7v + 5 < arg(A) < — n + — 5. This completes the proof. □ 

Theorem 4.5. Let m = 3 and let a e C m and < 5 < be fixed. Then 
W-i,i(a, A) =[-2uT* + O (\~ p )] exp [L(G\a), uT 2 A) - L{a, A)] 

- [2iw* + O (A~ p )] exp [-L{G 2 {a),u- 1 X) - L(a, A)] , 
as A — > oo in £/ie sector —5 < arg(A) < 7r — S. Also, 

W_i,i(a, A) =[-2iwf + O (A"")] exp [-L(a, uT 5 A) - L(tT 2 (a), uT 4 A)] 
+ [2wt + 0(A^)]exp [L(G(a),^ 3 A) -L(a,uT 5 A)] , 
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as A — > oo in the sector —n + 5 < arg(A) < 5. 

Proof. See Theorems 14 and 15 in [24] for a proof. □ 

5. Asymptotics of W-i, n (a, A) 

In this section, we will provide asymptotic expansions of W_i in (a, •), zeros of which will 
be closely related with the eigenvalues of H n . 
First, we treat the cases when 1 < n < |_ttJ. 

Theorem 5.1. Let 1 < n < [yj be an integer. Then W_i, n (a, •) admits the following 

asymptotic expansion 

(5.1) 

W_i, n (a,A) = [2cu 2 -^ L+ ^ Gn ~ 1 ^ +0(A-")]exp [L(G _1 (a), uj' 2 \) - L(G n -\a), w 2(n - 1} A)] , 
as A — >■ oo in the sector 

, 2(n — lW . _ , _ 4n7r 

(5.2) — + 5 < arg(A) < tt + 5. 

v; m + 2 - b\ j - m + 2 

Also, 

W_ 1>n (a,A) = [2^+^"^) + O (A"")] exp ^(G^a), uT 2 A) - L^^a), a; 2 ^- 1 ^)] 

(5.3) + [2w^ + ^ 1(a)) + O (\- p )} exp [L(G», w 2n A) - L(a, A)] , 
as A — >■ oo in the sector 

2(ra- 1)tt r ^ 2(n- 1)tt . 

5.4 i j! ^ < arg(A) < ^ -— + 5. 

y J m + 2 _ BU_ m + 2 

Proof. First we will prove (5.1) for the sector 

/r -\ 2(n — l)n . ^ ... ^ 4n7r 

(5.5) — + 5 < arg(A) < tt 5 

v 7 m + 2 -6V7- m + 2 

and the second part of the theorem by induction on n. 

The case when n — 1 is trivially satisfied by Theorem 4.3 and Corollary 4.4 since yu(a) + 

2v(a) = n{G-\a)). 

Suppose that (5.1) holds in the sector (5.5) for n — 1. From this induction hypothesis we 
have 

' 'H- : (G(- " 2 ' 

n-1 i ,,/nn-ll 



W 0) „(a, A) = w- 1 W_i,„_i(G(a),a; 2 A) 
(5.6) = [2^-^+^ G ^ 1(a)) + O (A"")] exp [L(a, A) - L^^a) V (n ~ 1} A)] , 



as A — )■ oo in the sector 



2(n-2)7r r ^ , 2u ^ 4(ra- 1)tt . 

+ 5 < arg(w 2 A) < tt ^ 5, 



m+2 ~~ ' — m+2 

that is, 

2n7r „ . . 4nn 

5.7 + 5 < arg A < tt - 5. 

v; m + 2 " 6W " m + 2 
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Also, from Lemma 4.1 if 1 < j < y + 1, then we have 

(5.8) W 0l , (a, A) = [2iaT* + O (A~ p )] exp [L(G'(a),u 2j - m ~ 2 X) + L(a, A)] , 
as A — > oo in the sector 

n — 1_ § < arg(A) < 7r — 5. 

m + 2 _ v ' ~ 

We solve (3.12) for W^_i(a, A) and set i = n to get 

/ KO x / u W-i,o(a, A)W„ lfc (a, A) W ,„(a, A)W_i )fc (a, A) 

(5.9) W_i,„(a, A) = — — — 1 — — — 

WV(a,A) Vy , fc (a,A) 

Set k = LfJ. Then since 1 < k - n < k = [f J , using (3.13), 

w , M _ 2^' 1 W)W ^ n (G'"(a),a; 2 "A) 1%> A)W , fc+1 (G'- 1 (a), a;' 2 A) 
1,nla ' J c«-iWo, fc (a,A) + w- 1 W 0>fc (a,A) 

2w MG- 1 W)[ 22a; - i ¥ i + 0(A^)]exp [L(G n (a),w 2fe - m - 2 A) + L(G n (a), u; 2 "A)] 

u;"- 1 ^-^ + O (A-")] exp [L(G fc (a), w 2fc -™- 2 A) + L(a, A)] 
[2a;- V+MG"- 1 ^)) + o ( A -p)] exp [ L ( a> A ) _ L(G n ~ 1 (a),u 2 ^ n ~ 1 ^\)] 

uj- l [2iuj-^ + O (A-")] exp [L(G k (a),u 2k - m ' 2 X) + L(a, A)] 
x [2iu~^ + O (A"")] exp [L(G k (a),u 2k - m - 2 \) + LfCT^a), uT 2 A)] 
= [2u-^ + ^ G ~ 1 ^ + (A"")] exp [L(G n (a),u 2n X)-L(a,X)] 

(5.10) + [2a;^ + ^ Gn_1 ^ + 0(A-")]exp [L(G _1 (a), uT 2 A) - L(G n -\a), u 2 ^ X)} , 

where we used (5.6) for Wo, n (o, A) and (5.8) for everything else, provided that A lies in (5.7) 
and that 

_ 4 (L^J - rx) tt < (w2 „ A) < _ 
m + 2 _ v 7_ 

41 - I-7T 

7T - + 5 < arg(A) < 7T - 5 

m + 2 - _ 

^ _ 4 ( If J + 1) tt (ar 2 A) < ^ _ s 



that is, 



2n7r 4n7r 

+ 5 < arg(A) < tt 5. 



m + 2 ' m + 2 

Thus, the second part of the theorem is proved by induction. 
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Next in order to prove the first part of the theorem for the sector (5.5), we will determine 
which term in (5.10) dominates as A — > oo. To do that, we look at 

Re (L(G- 1 (a),w- 2 A) -L(G n - 1 (a),w 2(n ~ 1) A)) - Re (L(G n {a), u 2n \) - L{a, A)) 



= K n 



, 2tt m + 2 

cos I 1 arg( A) — cos 

m 2m 

2mv m + 2 



2(n-l)n , m + 2 

+ — ; arg(A) 



m 



2m 



— I cos I 1 — arg(A) — cos 



2K m sin 



m 2m 



m + 2 
2m 



arg(A) 



| A | 3+™(l + (l)) 



(5.11) = 4iir m siii 



m 



n-zr 



m 



sm 



>-2)tt m + 2 

+ ^T arg(A) 

n7r m + 2 



m 



+ sin h 



7T 

cos ( — ) sm 

m. 



m 2m 
(n — l)7r m + 2 



arg(A) 



|A|*+™(1 + (1)) 



m 



+ ^ r arg(A))|A|^(l + o(l)), 



that tends to positive infinity as A — > oo (and hence the second term in (5.10) dominates) if 



2(n-l)7r 



+ 5 < arg(A) < ir 



4nn 



5. 



m+2 1 — "^bV'V — m +2 

We still need to prove (5.1) for the sector 



4n7r 4mr 

7r r — < arg A < n h 0, 

~ v ' ~ m + 2 



(5.12) 

v ; m + 2 

for which we use induction on n again. 
When n = 1, (5.1) holds by Lemma 4.4. 

Suppose that (5.1) in the sector (5.12) for n - 1 with 2 < n < [f J. Then (3.13) and (5.9) 
with k = n + 1 yield 

1/W , n 2u^ G - 1 ^W 0A (G n (a),u 2n X) W- 1 , n - 1 (G(a),o; 2 A)Wo, ra+2 (G- 1 (a)^- 2 A) 



w n " 1 >Vo ) n+i(a, A) 



W 0) „+i(a, A) 



If < arg(A) < ^-^ + 5, themr-^^-S < arg(w 2 A) < tt-%^ + 5. So 



m+2 



(5.13) 

v ( G -i(a)) + ^(.)) W_ 1 , n _ 1 (G(a),a; 2 A)Wo, ra+2 (G- 1 (a),a;- 2 A) 
linla ' J w»- 1 Wo > ^i(a,A) + W 0>B+ i(a,A) 

Au v(G-Ha))+v(G"(a)) 



+ 



u"- 1 [2ioj-* + O (A-")] exp [L(G n+1 (a) , u 2n ~ m X) + L(a, A)] 
[2 W 3?+MG»- 1 W) + ( A -p)] exp [ L ( a> A ) _ L(G n ~ 1 (a), w 2(n_1) A)] 



[2?cj-^ + O (A-")] exp [L(G n+1 (a),w 2n - m A) + L(a, A)] 
x pic;" 2 * 2 + O (A"')] exp [L(G n+1 (a),a; 2n - m A) + L(G- 1 (a),^- 2 A)] 
= [_2^^ + ^ G ~ 1 ( a )) + ^ G "( a » + O (A"')] exp [-L(G n+1 (a)V n - m A) - L(a, A)] 
+ p^+MG"- 1 ^)) + ( A -p)] exp [ L ( G -i( a ) }W -2 A ) _ L(G"- 1 (a),a; 2 ("- 1 )A)] , 
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where we use the induction hypothesis for W_i >n _i((j(a), a; 2 A) and use (4.1) for Wo, n +i( a , A) 
and Wo, Tt +2(G'" 1 (a), u~ 2 \). Next, we use an argument similar (5.11) to complete the induc- 
tion step. Thus, the theorem is proved. □ 

Next we investigate Wo,\_m.\(a, A). 

Theorem 5.2. If m> A is an even integer, then 
W_i, LfJ (a, A) 

(5.14) = _[2 CJ 2+MG- 1 W)+MG^(a)) + ( X - P y exp [_ L ( G f + 1 ( a ), A) - L(a, A)] 

m — 2 r _„ 

- [ 2w 2 +M«)+MG^(«)) + o (x-p)\ exp -L(G m ^{a),u m - 2 X) - L{G m {a) , u m ~ 2 X) 
as A — >■ oo in £/ie sector 

Att An 

(5.15) -ttH <5 < arg(A) < -7T H h 5. 

v; m + 2 - &w - m + 2 

If m > A is an odd integer, then 

W_ UfJ (a,A) = [2u4 + 0(A-")]exp [-L(G^(a), a; -1 A) - L(a, A) 

(5.16) + [2a;f + 0(A- p )]exp \L(G m+ \a),uj- 2 X) - L(G^ {a),u m ^X) 
as A — >■ oo in the sector 

4 71 " r „ , n , 67T 

-7T + — + S < arg A < -7T + — + 5. 

m + 2 m + 2 

Proof. We will use (5.13) with n = [f J, that is, 

4a; MG- 1 («))+MG l ^ J (a)) W_ 1)Lf J _ 1 (G'(a),a; 2 A)>Vo, Lf J+2 (G" 1 (a), w" 2 A) 
W_Uf J A) = uj^W Ufi+1 (a,X) + ~ Wo, Lf j +1 (a,A) ' 

When m is even, say m = 2k, 

Wo,k + 2(G-\a),uj- 2 X) = W m+ 2,k + 2(G-\a),uj- 2 X) 

= -u-^W-^iG^ia),^^) 
= u- 2 W- l ^ 1 {G k+2 {a),u 2 \). 

So 

^(g-^a^+Mg^a)) W_ 1 , fc _ 1 (C(a),a; 2 A)Wo, fc+2 (G- 1 (a),a;- 2 A) 
1>fcla ' } u^W , k+1 (a,X) + W 0ife+1 (a,A) 

_ ^(G-^+MG^a)) W-i, fe -i(G(a), a; 2 A)W-i, fc -i(G fe + 2 (a), a; 2 A) 
> 'VV,,./... ; U/.A) ^-'VV,,./.-. :l".A) 

Since A lies in (5.15), 

2 a?j- 2 )' +f <- ir+ -gL ; - f < ttg(M . A) < ir - 4 (Lgj-i)T +f , 



m + 2 _ m + 2 _ ov y_ m + 2 
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yV-!, k (a, A) 

j^^Ml+M^W) ^ W- 1 , fc - 1 (G(a),a; 2 A)Wo, fc+2 (G- 1 (a)^- 2 A) 
w*- 1 >Vo ) fc+i(a,A) ~ ~ W ,fe+i(a,A) 
Au ^G-Ha))+^GHa)) W-i |ife -i (g(a) , i^ 2 A) W-i, fc _i (G fc+2 (a) , CO 2 A) 
w fc - 1 W ,fe+i(a,A) + w 2 W , fc +i(a,A) 

4w MG- 1 (a))+^(G fc (a)) 

~ w*- 1 ^-** 1 + O (A-")] exp [L(G fe+1 (a),a; 2fe - m A) + L(a, A)] 
pc^+M^-V)) + o (A"')] exp [L(a, A) - L(G fe - 1 (a), co 2 ^X)] 
u 2 [2iu-^ + O (\-p)] exp [L(G k + l {a), u 2k ~ m X) + L(a, A)] 
x [2^+MG 2fc M) + ( A -„)j exp [ L (G fc+1 (a), A) - L(G 2k (a), u 2 ^ A)] 

= [_2 ?w ^ + ^ G "( a )) + ^ Gfc ( a )) + O (A"")] exp [-L(G fc+1 (a), A) - L(a, A)] 

- p^+MG^WH^W) + O (A"')] exp [-L(G fe - 1 (a),a; 2 ( fe - 1 )A) - L(G 2k (a) , u 2 ^ X)} 
= _[2 ?w ^+MG- 1 (a))+M(G'=(a)) + q ^] exp ( a ) , A) - L(a, A)] 

- p^^+MW+M^-W) + o (\-p)\ exp [-L(G fc - 1 (a),a; 2 ( fc - 1 )A) - L(G 2k (a) , u 2(k ^ X)} , 

where we used (5.1) for W-i tk -i(G(a), •) and W-i tk -i(G k+2 (a), •), and (4.1) for W 0)fc+1 (a, •). 
Finally, we use a;" 21 ^ = — i, to get the desired asymptotic expansion of W-i ! [m.\(a : A). 
Next we investigate the case when m is odd, say m = 2k + 1 (so = k). 

W 0)Lf j+i(a, A) = a;- 1 !^-!,^ j (G(a), w 2 A) 

and 

W , fc+2 (G- 1 (a),a;- 2 A) = W^^G^ta), uT 2 A) 

= -a;- fe - 2 Wo, fc+1 (G fe+1 (a),a; 2fc+2 A) 
= a;^W , fc+1 (G' fe+1 (a),a;- 1 A). 

Similarly to the proof of the theorem for m even, 

4 ^->)) + MG fc (a)) W_ 1 , fc _ 1 (G(a)^ 2 A)>Vo, fc+2 (G- 1 (a),a;^A) 
lMfl ' j w*- 1 W 0>fc+ i(a,A) + Wo, fc+1 (a,A) 

_ ^(G-Ha)) + ,(GHa)) ^ yy_ lk _ l (G(a),u 2 X)W 0ik+ i(G k+1 (a),uj' 1 X) 
w fc - 1 W 0l fc+i(o,A) w^W ,fc+i(a,A) 
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4w M(G- 1 («))+^(G fc (a)) 



Lu^ftiuj-^ + O (X-p)} exp [LiC+^a), u 2k ~ m X) + L{a, A)] 
^sjL+rtak-Ha)) + ( A - P )] exp [ L ( a> A) _ L (G k -\a), u 2 ^X)] 

u^ftiu"^ + O (X~p)] exp [L(G k+1 (a),u 2k - m X) + L(a, A)] 

x [2iuT^ + O (A"")] exp [L(G 2k+2 (a), oo^'^X) + L(G k+1 (a) , u' 1 X)} 

= l^icu^^- 1 ^^^ + O (X-p)] exp [-L(G k+ \a), cu^X) - L(a, A)] 

+ [2u 2 ^+^ Gk - 1 ^ +C(A-")]exp [L(G 2fc+2 (a),^- 2 A) - L(G fe - 1 (a), w 2 **- 1 ^)] 

= -[2iuj^ + f + O (A"')] exp [-L(G k+1 (a), co^X) - L(a, A)] 

+ [2o;^ + t + O (A"")] exp [L(G 2k+2 (a) , u~ 2 X) - L{G k -\a), u 2 ^X)] , 

where we use (5.1) for W_i,fc_i((3 f (a), •), and use (4.1) for Wo,jt+i(a, •) and Wo > k+i(G k+1 (a), •). 
Finally, we use w 12 ^ 2 = z, to get the asymptotic expansion of W-i^sj (a, A). This completes 
the proof. □ 

The order o/ an entire function g is defined by 

log log M(r, a) 
lim sup : , 

r^oo log r 

where M(r,g) = max{|o(re j61 ) | : < 8 < 2n} for r > 0. If for some positive real numbers 
a, ci, C2, we have exp[cir CT ] < M(r,g) < exp[c2r CT ] for all large r, then the order of g is a. 

Corollary 5.3. Let 1 < n < y. Tnen tne entire functions W-i, n (a, •) are o/ order § + ^, 

and aence tney nave infinitely many zeros in the complex plane. Moreover, W-i, n (a, •) nave 

at most finitely many zeros in the sector 

2(n - 1W r ^ ... 47r 

+ 5 < arg(A) < tt H 5. 



m + 2 ~ ~ m + 2 

6. Relation between eigenvalues of H t and zeros of W_i j7l (a, •) 

In this section, we will relate the eigenvalues of Hi with zeros of some entire function 
W-i,„(a,.). 

Suppose that £ = 2s — 1 is odd with 1 < £ < m — 1. Then (1.1) becomes (3.2) by the 
scaling v(z) = u(—iz), and v decays in the Stokes sectors S- s and S s . Since / s _i and / s are 
linearly independent, for some D s and D s one can write 

f- s (z, a, A) = D s (a, A)/ a _i(z, a, A) + 5 s (a, A)/ a (z, a, A). 

Then one finds 

LUa, A) = — — —. — rr and DJa, A) = — — - — -, — -— . 

V ; W s _ M (a,A) sV ; W S;S _i(a,A) 

Also it is easy to see that A is an eigenvalue of Hi if and only if D s (a,X) = if and only 
if W_ s , s (a,A) = 0. Since W_ SjS (a,A) = u a - 1 W- 1 , 2a -i(G- a+1 (a), u~ 2s+2 X), by Corollary 5.3 
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W- S)S (a, A) has at most finitely many zeros in the sector — 2 ^^ n + ^ — ar g( w 
71 + ^+2 ~ ^ that is ' 

/n , 2(2s + 1)tt . 

5 < arg(A) < 7T + — 5. 

~ v ' ~ m + 2 

Next, by symmetry one can show that W- SjS (a, A) has at most finitely many zeros in the 
sector 7r < arg(A) < 2n — 6. For that, we examine H e with P(z) replaced by P(z) whose 
coefficient vector is a := (a l5 a 2 , . . . ,a m ). Then one sees that W_ SiS (a, A) = if and only if 
W_ s>s (a, A) = 0. W- S)S (a, A) has at most finitely many zeros in the sector 5 < arg(A) < n 
by the arguments above. Thus, W_ SjS (a, A) has at most finitely many zeros in the sector 
7r < arg(A) < 2n — 5, and has infinitely many zeros in the sector | arg(A)| < 5 since it is an 
entire function of order | + ^e(0,l). 

Suppose that £ = 2s is even with 1 < £ < m — 1. Then (1.1) becomes (3.4) by the scaling 
y(z) = u(—iu~^z), and y decays in the Stokes sectors S- 8 and S s+1 . We then see that the 
coefficient vector a of the polynomial uj~ 1 P{uj~^z) becomes 

a = G~z(a). 

Now one can express f- s as a linear combination of f s and f s+ \ as follows. 

fs(z,a,u A) = — fs{z,a,u A) + — — / s+1 (^,a,o; A). 

yV s , s+ i(a,u; U) >V s+ i iS (a,cj 

Thus, A is an eigenvalue of iJ^ if and only if W- s , s +i(a, u; _1 A) = 0. Since 

W-^+ifoa;-^) = a; s_1 W_i ) 2 S (G' _s+1 (a),a; _2s+1 A), 

by Corollary 5.3, W_ SiS+ i(a, u; _1 A) has at most finitely many zeros in the sector - 2(2 ^ )7r + 
5 < arg(cj- 2s+1 A) < tt + ^ - o", that is, 

5 < arg(A) < tt H 5. 

~ v 7 _ m + 2 

This is true for each a G C m . So one can show that W_ S)S+ i(a, u; _1 A) has at most finitely 
many zeros in the sector n < arg(A) < 2-k — S by symmetry, similar to the case when £ is 
odd. Thus, W-a jS +i(o, A) has infinitely many zeros in the sector | arg(A)| < 8. 

7. Proof of Theorem 1.1 when 1 < £ < 

In this section, we prove Theorem 1.1 for 1 < I < [yj and in doing so, we will use the 
following proposition on univalent functions. 

Proposition 7.1 ([17, p. 216]). Let A{jj) be analytic in the region S = {fi G C : a < 
arg(/x) < j3, \fx\ > M} for some a, f3 G M. with f3 — a < n and for some M > 0. Suppose that 
A(/j) — > as n — > oo in S. Then there exist a, {3, M such that /x(l + A(fj)) is univalent in 
{/iGC:5< arg(/i) < f3, > M} C S. 

Proof. See [17, p. 216] or [27, Thm. 3.4] for a proof. □ 
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We will consider two cases; when £ is odd and when £ is even. 

Proof of Theorem 1.1 when 1 < I < [yj is odd. We will closely follow his proof of Theorem 
29.1 in [27] where Sibuya computed the leading term in the asymptotics (1.3) for £ = 1. 

Suppose that 1 < £ = 2s — 1 < [yj is odd. Recall that when £ is odd, A is an eigenvalue 
of Hi if and only if W_ S)S (a, A) = 0. Also, in Section 6 we showed that W_ SjS (a,A) has all 
zeros except finitely many in the sector | arg(A)| < 5 and hence, all the eigenvalues A of Hi 
lie in the sector | arg(A)| < 5 if |A| is large enough. 

Since 



W- S;S {a, A) = u^W-^s-^G'^ia), 



-2s+2 



A), 



we will use (5.3) to investigate asymptotics of large eigenvalues. Suppose that W- SjS (a, A) = 
and |A| is large enough. Then from (5.3) with n — 2s — 1, and with a and A replaced by 
G~ s+1 (a) and u~ 2s+2 \, respectively, we have 



Anv{G s {a)) 
m+2 



[1 + (\- p )] exp [L(G s (a),u 2s X) - L(G' s (a), uT 2s A)] 

x exp [L(G s -\a),u 2s - 2 X) - L(G~ s+ \a),oj- 2s+2 \)} 

Also, since [1 + (A"")] = exp [O (A"')] and u,MG») = exp 

exp [L(G s (a),u 2s X) - L(G~ s (a), w~ 2s A)] 

(7.1) x exp 

For each odd integer £ = 2s — 1 in 1 < £ < [_ttJ , we define 
h m/ {X) =L(G s (a),u 2s \) - L(G~ s (a),u~ 2s \) 



L(G s ^ 1 (a),o; 2s " 2 A) - L(G- s+1 (a), u 



2s+2 _ 4*v(G'(a)) Q (A _ P) 
m + 2 v ' 



= -1. 



+ L(G s -\a),co 2s - 2 X) - L(G- s+l (a),u- 2s+2 \) - 4 ™} G ° < ? )) + O (A"') . 



Then by Corollary 3.4, 



= 2iK, 



1st: , 



m i e ™ -e 

2S7T 

sin 

m 



2an ; 2(3-1)71 
m — |— g m 



2(3-1)71 



— e 



m + 2 



\— (l + o(l)) 



+ sin 



2sir 



m 



2TT \ \ , m+2 . 

A" l + o 1 

m ' ' 



(7.2) 



a-ts f 7l \ ■ A 2 « - l)7r\ x 12+2 

= 4zif m cos — sin — A 2m (1 + o(l)) as A ->■ 00, 

Vm/ V ml 



in the sector | arg(A)| < 5. Since K rn > and < ^ 2s ~ 1 ' 7T < 7^ the function h m ^{-) maps the 
region |A| > Mi for some large Mi and | arg(A)| < 5 into a region containing |A| > M 2 for 
some large M 2 and | arg(A) — || < E\ for some £1 > 0. Also, h m ^{-) is analytic in this region 
due to the analyticity of f(0, a, A) in Theorem 3.2. Following Sibuya, we will show that for 
every large positive integer n there exists A n such that 

(7.3) h m/ {X n ) = (2n + 1) m. 
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Next, Proposition 7.1 with /i = A 2 ^ 2 implies that there exist M{ > and < 5' < | such 
that fj,(l + A(fj,)) is univalent in the sector | arg(yu)| < 5' and > M[. Thus, for each n£Z, 
there is at most one fi n in the sector | arg(yu)| < 5' and > M{ such that 

(7.4) // n (l + A(/x n )) = B n , 

where A is the error term in (7.2) and 

(2n + 1)tt 



-Br; 



4K m cos (^) sin 



(2s-1)tt 



Then there exists M" > M[ such that in the sector | arg(/i)| < 5' and > M", 

\m\<\- 

Also we can take n large enough so that \B n > M". Next, we define 

A n :={\A^)\:\wg{n)\<5',\n\>B n }. 

Then lim^oo A n = 

Suppose that n is so large that A n <\ sin(<5') < |. Then the disk defined by \/i — B n \ < 
i-2A is contained in the sector | arg(yu) | < 5' and > \B n . Moreover, on the circle 
I u — B I = A " B 

\h L ^nl 1-2A„ n > 

\n - B n \ - \fi\\A(n)\ > \n - B n \ - \n\An 

> \n - B n \(l - An) - A n B n 

> \Z 2 A AnBn ~ AA > °' 

Thus, by the Rouche's theorem in complex analysis (see, e. g., [7, p. 125]), fi — B n and 
fi(l + A(n)) — B n have the same number of zeros in the disk and hence, (7.4) has exactly 
one fj, n in the sector. Therefore, there exists iV = N(m, a) > such that (7.3) has exactly 
one solution A„ for all integers n > N. 
Next, since 

— — ln(aT s A) + — — Info; 2s A) 

m m 

HG-Ha)) ^ HG-'(a)) ln( ^ x) _ MO)),, 
m v ' m v ' m + 2 



m m + 2 

= - —u(G s (a))ni, 
m 

using Corollary 3.3, (7.3) becomes 

m+l . . 

(7.5) (2n + 1) ni = e^{a)\l + ^ v(G s (a))ni + O (\~ p ) , 

j=0 
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where 

e,» =X m ,^(a))u; 2s (H^) -K mj (G'- s (a))a;- 2s (5+^) 

+ K m ,(a) (w-^+^KH^) _ w (-iW-2(.-D(H^)) 

where we used Lemma 3.5. 

In summary, we have showed that for each a G C m , there exists iVo = N (m, a) < N such 
that {A n }„>jv is the set of all eigenvalues that satisfy (7.5), since there are only at most 
finitely many eigenvalues outside | arg(A)| < e or inside |A| < M for any s > 0, M > 0. In 
order to complete proof, we need to show that N does not depend on a G C m . 

We will prove that for each R > 0, if \a\ < R, then N (m, a) = N (m, 0). The error terms 
in (3.5) and (3.6) are uniform on the closed ball \a\ < R and hence, so is A(y) as fx — > oo in 
the sector | arg(//)| < 5'. Thus, we can choose independent of \a\ < R. Since //(l + A(/j,)) 
is univalent in the sector | arg(/i)| < 5' and > M" > M[, by the implicit function theorem 
and (7.4) the function a >->■ /i n (a) is continuous on |a| < i? for each n > N and hence, so is 
a i — y A„(a) for n > N. 

Next, we claim that for all |a| < i?, there are exactly the same number of eigenvalues that 
are not in {\ n (a)} n >N- To prove this, we will use the Hurwitz's theorem (see, e. g., [7, p. 
152]) in complex analysis. That is, if a sequence of analytic functions converges uniformly 
to an analytic function on any compact sets, then eventually functions in the sequence and 
the limit function have the same number of zeros in any open set whose boundary does 
not contain any zeros of the limit function. The Hurwitz's theorem implies that since the 
eigenvalues are the zeros of the entire function W_ SjS (a, A), they vary continuously as a 
and hence, there is no sudden appearance or disappearance of eigenvalues. Also, none of 
the eigenvalues that are not in {X n (a)} n >N can be continuously mapped to A n (a) for some 
n > N as a varies and hence, the claim is proved. This completes the proof. □ 

Next we prove Theorem 1.1 for 1 < i < |_ttJ is even. 

Proof of Theorem 1.1 when 1 < t < [yj is even. Proof is very similar to the case when 1 < 
£ < [f J is odd. Let £ = 2s for some s G N. 



23 

Recall that A is an eigenvalue of H e if and only if W_i 5 2 S (G~ s+1 (a), u~ 2s+1 X) = 0. Then 
from (5.3), we have 

exp [L(G s+1 (a),co 2s+1 X) - L(G~ s (a), oT^A)] 

(7.6) x exp [L(G 8 (a),u} 2s ~ 1 X) - L((r s+1 (5), uT 2s+1 A) + O (A"')] = -1, 

where a = G~^(a) and where we used [1 + (A~ p )] = exp[0 (A~ p )] again. 
Like in the case when 1 < £ < [yj is odd, from Lemma 3.3, we have 

m+i x ^ 

(7.7) (2n + 1) m = ^ e^(a)Al + ^ + O (A"') , 

j=o 

where for < j < m + 1 , 

etj(a) =K m ,{G s+ Ha))uj^ +1 ^ +1 -^) - /^-(G— i(o)) W -^ +1 >(H^) 

+ ^((Ti^H^) - X m ,(G'- s ^(a))a;-( 2s - 1 )(H^) 

= £(-l)^A*(«) 

fc=0 

=4^(-l) fc K m ,, fe ^(a)sin ffi^QM cos (2-^) , 
k=o \ m j \ m j 

where we used Lemma 3.5 and £ = 2s. Then we use the arguments in the case when £ is odd 
to show that there exists N = N (m) such that the set {A n } n >jv of all eigenvalues satisfy 
(7.7) and hence, the proof is completed. □ 

8. Proof of Theorem 1.1 when £ = \J^\ and when f < £ < m - 1. 

In this section, we prove Theorem 1.1 for £ = [_ttJ • We first prove the theorem when m is 
even, and later, we will treat the cases when m is odd. Then at the end of the section, we 
will prove the theorem when y < £ < m — 1, by scaling. Proof of existence of iVo = N (m) 
is the same as in Section 7, so below we will omit this part of proof. 

8.1. When m is even. We further divide the case into when £ is odd and when £ is even. 

Proof of Theorem 1.1 when m is even and £ = [yj is odd. Suppose that m is even and £ = 
Y = 2s — 1 for some s£N. Recall that A is an eigenvalue of Hi if and only if W_ SjS (a, A) = 
if and only if 14/-i, 2s -i(G- s+1 (a),w- 2s+2 A) = 0. 

So if W- ll 2s-i(G- 8+1 (a),u- 28+2 \) = and if |A| is large enough, then from (5.14), 

exp [L(G 3s ^(a),u 2s - 2 X) - L(G S+1 (a) , u' 2s+2 X)] 

x exp [L(G s -\a),u 2s ~ 2 X) - L(G- s+1 (a),cu~ 2s+2 X) + O (A"")] = -^ G >». 
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Then by Lemma 3.3, 
f8v(G s (a)) 



V m + 2 



m+l 



j=0 



m(m + 2) 



where for < j < m + 1, the coefficients eej(a) are given by 

etj{a) =K mj {G^-\a))^ s - l ^ + ^) - K mJ (G s+ \a))u;- 2 ^(^) 

+ K mJ (G s -\a))^ s -^ +1 -^) - K m j{G- s +\a))u-^-^^) 
=K mJ (a) ( w MW+2M)(i+bi) _ u -(^-2M)(i+3 S i) 

+a; - (s -l) j+ 2( S -l)(I + ^) _ w ( s _i)j_ 2 ( s _i)(i + ^ 

'(1-^(23-2)^ 



--2iK mJ (a) (l + (-l) J ) sin 



m 



AiK m Msm[^-J^] cos^ (1 - j)7r 



where we used Lemma 3.5. 



□ 



Proof of Theorem 1.1 when m is even and £ = [ttJ is even. Let £ = y = 2s for some s G N. 
Recall that A is an eigenvalue of H e if and only if W-i : 2 S (G~ s+1 (a), uj~ 2s+1 \) = 0. 

Suppose that W_i, 2s (G'" s+1 (a), u~ 2s+1 X) = W_i, 2s (G^(a), ^- 2s+1 A) = 0. Then by 
(5.14), 



exp 
x exp 

Then like before, 



L(G 3s+ ^(a),a; 2s - 1 A) - L((G s+ § (a), u' 2s+1 X) 

L(G s ^^( y a),u 2s ^ 1 X) - L(G~ S+ ^ (a) , u~ 2s+1 \) + O (A~ p ) 



= -1. 



m+l 



(2ra + 1) wi = J2 e£,i(a)Al +1 ™ + O (A"') , 
i=o 

where for < j < m + 1 , 

+ K^G'-^a))^-^^) - K m ,{G-°+Ha))u-^- l ^+^) 



=2iX)(-l)*A mj - >fc 6 J -fc(o) (! + sin 



fc=0 



(1-J)(2S-1)7T 



m 



\>Yj l)'/^,./,/',/.!")^.. ( cos ^ • /) " 



fc=0 



□ 
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8.2. When m is odd. We divide the case into when £ is odd and when £ is even. 



Proof of Theorem 1.1 when m and £ = [yj are odd. Let m and £ = rn ^- = 2s — 1 be odd. 
Suppose that A is an eigenvalue of Hi. Since A is an eigenvalue of Hi if and only if 
W-i,2 S -i(G~ s+1 (a), u~ 2s+2 X) = 0, if |A| is large enough, then from (5.16), 

exp [L{G s -\a),u 2s - 2 \) - L(G- s+1 (a),uj- 2s+2 X)] 

x exp [-L(G 3s+1 (a),uj- 2s \) - L(G s+ \a), uT 2s+1 A) + O (A"')] = -1. 

Then, from Lemma 3.3, 

m+l 



(2n +l)m = J2 eij(a)\l +1 ™ + O (X~ p ) , 

3=0 

where for < j < m + 1 , 

e,» =K m ,{G s -\a))^ 2s - 2 ^^) - K m ,{G- s+ \a))u^ 2s+2 ^ + ^) 
- K rn ,{G^\a))uj- 2 <^) - K m ,{G s+ \a))J- 2s+l ^ +1 ^) 
=K m ,{a) (u-^ s - l ^ 2s - 2 ^ +1 ^) - u^- {2s - 2 ^ +1 -^) 

_ u; -3(3s+l)-2s(l + ^) _ w - j(s+ l)-(2 s -l)(I + i^)^ 

=K md (a) (u^ i2s - 2 ^ + ^) - u^ 1 -^- 2 ^) 

n . T , t \ ( • f {m - 3)(1 - j)ir\ . /(m+l)(l-j> 
--2iK mJ (a) ( sin ( ^ — ) + sin ' 



2m j \ 2m 



-4iK mi j(a) sin ( | cos 



ml \ m 



7T 



□ 



Proof of Theorem 1.1 when m is odd and £ = [^J even. Let £ = = 2s for some s e 
N. Then A is an eigenvalue of # £ if and only if W_i, 2s (G" s+1 (a), uj- 2s+1 X) = 0. 
If W_i, 2s (G- s+1 (a),cj- 2s+1 A) = and if |A| is large enough, then from (5.16), 

exp [LiG^a)^ 23 -^) - L(G- s+1 (a),uj- 2s+1 X)] 

x exp [-L(G 3s+3 (a),u' 2s ^ 1 X) - L(G s+2 (a), uT 2s A) + O (A"")] = -1. 
Then, from Lemma 3.3, 

TO+l 



(2n +l)ni = J2 e^(a)Al +1 ™ + O (A"") , 

3=0 
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where for < j < m + 1 , 

e,» =^-(G-i(fl)) W M^) - ^((^(aKM^) 
- K mj (G 3s+ f (a))c-( 2s+1 Ki+^) - A^(G*+I(a))uT 2s (^) 

=2* 2J-l)*tf mJl A fc (a) sm ( — — 1 + sm ( — — 

fc=o L ^ 

=4i£(-l)'tf m , M W«) sin cos 



m / \ m 



□ 



Theorem 1.1 for 1 < £ < y has been proved. Next we prove Theorem 1.1 for y < £ < 
m — 1, by the change of the variables z i— >■ — z. 

8.3. When f <£<m-l. 

Proof of Theorem 1.1 when y < i < m — 1. Suppose that £ > y. If w is an eigenfunction of 
Hi, then t>(z) = u(—z, A) solves 

+ [(-l)- € (-iz) m - P(-iz)] v(z) = Xv(z), 

and 

/x ii ii | ((m — £) + 1)7T 

v(z) — » exponentially, as z — > oo along the two rays argz = — — ± . 

The coefficient vector of P{—z) is ((— l)" 1-1 ^, (— l) m ~ 2 a 2 , . . . , (— l) 1 a m _i, a m ). Certainly, 

. f (l-j)( m -£)n \ . f(l-j)£ir 

sm = (— l) J sm ' 



m J \ m 

Also, one can find from Lemma 3.5 that for < k < j, 

6j- >fc ((-l) m-1 ai, (-l) m_2 a 2 , . . • , -a m _i, a m ) = (-l) mk ^b jjk (a u a 2 , . . . , a m _i, a m ). 

Moreover, c m -e !j ((-i) m ~ 1 a 1 , (-l) m ~ 2 a 2 , . . . , -a m _i, a m ) = Qj(ai, a 2 , . . . , a m _i, a m ). This com- 
pletes proof of Theorem 1.1. □ 
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